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Abstract 

We diagonalise the transfer matrix of boundary ABF models using bosonized vertex opera- 
tors. We compute the boundary S-matrix and check the scahng hmit against known resuhs for 
perturbed boundary conformal field theories. 



CJ 
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■ 1 Introduction 



The theory of solvable lattice models, which originated in the work of Bethe and Onsager in the 30- 
40's, has matured in the last 20 years. Various models were found to be solvable and various methods 



were invented to solve these models. In the course of this work, some remarkable interactions 

O . 

. occured between the theory of solvable lattice models and other branches of mathematics and 



physics, e.g., representation theory and conformal field theory. 
0^ . The ABF models we are going to study in this paper were introduced in Baxter's book [Q] 

presented the corner transfer matrix method, by which the one point correlation functions were 
computed. The general results obtained in ||l| strongly suggested a link between solvable lattice 
^ ■ models and the representation theory of the Virasoro algebra. 

Such a link was manifest from the very beginning in conformal field theory. The representation 
' theory of the Virasoro algebra was very successfully applied to conformal field theory to obtain 
^ ■ correlation functions and the differential equations which they obeyed. 

^ Solvability of lattice models was understood by means of the method of commuting transfer 

matrices. The algebraic structure behind it was given shape by the idea of g-deformation. Finally, 
in the corner transfer matrix method was given its correct place in the representation the- 

ory. The link was established between the highest weight representations of the quantum affine 
algebra Uq{sl2) and the XXZ model in the anti- ferromagnetic regime. Namely, the spaces of the 
eigenvectors of the corner transfer matrix of the latter were identified with the level 1 highest 
weight representations of the former, and the half transfer matrices were identified with the inter- 
twiners of these representations. In fact, a similar structure had been known in conformal field 
theory. The operators of this kind were called vertex operators. These operators were realized 
using representation theory [0. The method of bosonization was also very powerful P]. 
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As for the ABF models, such a scheme was only recently used successfully . The bosonization 
method of conformal field theory, by which the minimal unitary models were solved, was deformed 
to give a realization of the half transfer matrices of the ABF models. 

The boundary problem has been studied in both lattice and continuum theories. In [Q, the 
vertex operator method was extended to the boundary XXZ model. The aim of this paper is to 
extend the bosonization method to the boundary ABF models. 

We start from a known solution of the boundary Yang-Baxter equation . We then introduce 
the boundary transfer matrix of Sklyanin-type in terms of a product of vertex operators. Our goal 
is to diagonalise the boundary transfer matrix. We shall give an explicit formula for the boundary 
vacuum state and compute the boundary 5-matrix. The bosonization scheme of 0] is based on 
a cohomological construction. Thus, in the boundary ABF models, the physical spaces of states 
are realized as sub-quotients of the bosonic Fock spaces on which the screening operators act. 
One of the main points in our work is to show that the boundary vacuum states belong to the 
sub-quotients. 

The scaling limit of the bulk ABF models in regime III is described by the (p^i^s) deformations 



of the c = 1 — 6/r(r — 1) rational conformal field theories. In |10] and [^T|, boundary S-matrices 
are computed for the perturbed r = 4 and r = 5 models on a half line (the Ising and Tri-critical 
Ising models respectively). We check that our boundary S-matrices agree with these results in the 
scaling limit. 

The plan of the paper is as follows. Section 2 prepares the boundary Boltzmann weights and 
the transfer matrix. Section 3 is a summary of the bosonization scheme. We give the boundary 
vacuum state in Section 4, and compute the boundary 5-matrices in Section 5. Section 6 takes the 
scaling limit. Appendix A gives useful formulas for operator products. Appendix B gives a proof 
of the eigenvector relation for the boundary vacuum states. 



2 Bulk and boundary weights 

We recall the Boltzmann weights of the ABF model, and set up the boundary transfer matrix in 
terms of the vertex operators. 

2.1 The bulk weights 

The ABF model has two parameters x and r. We assume that < x < 1 and r ^ 4 (r G Z). We 
use the symbol [u] for the elliptic theta function. 

Oqiz) = {z;q)oo{qz'^;q)oo{q;q)oo, 

OO 

{z;qi,- ■ ■ ,Qm)oo = n (1 - • • -^m 
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There are six kinds of configurations around a face and the corresponding Boltzmann weights are 
given as fohows Q. 



W 



W 



w 



k±l 


k 


k 




k 


k±l 


k±l 


k 


k 


k±l 




k 



u = R{u) 



u = R{u] 



u = Riu 



[i-u][ky 

[k^l\[u] 



[k][i-uy 

Here k is an integer such that 1 ^ A; ^ r — 1, and the normahsation factor R{u) is given by 
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The Boltzmann weights satisfy the following relations. 
• Yang-Baxter equation 
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2.2 The vertex operators 

Consider the corner transfer matrices A(Q,B(C,), C{C), D{(^) which represent the NW, SW, SE, NE 
quadrants, respectively. Let Ci^k be the space of the eigenvectors of A{C,) in the sector such that 
the central height is equal to k and the boundary heights are (/,/ + 1). We denote by $('^+'^'^)(^) 
the half-infinite transfer matrix extending to infinity in the north. This is an operator 

Similarly, we denote by the half-infinite transfer matrix extending to infinity in the 

west. This is an operator 

They satisfy the following relations. 
• Exchange relation 
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b c 
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Duality 



• Inversion relation 
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(2.9) 



These operators are realized as vertex operators acting on the bosonic Fock space J^i^k in Section 
3.1. 



2.3 The boundary weights. 



We follow Sklyanin's scheme in dealing with the boundary ABF model |12]. Boundary weights 
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are given to the boundary configurations. We restricts to the diagonal case, e = e' , and use a 
solution of the boundary Yang-Baxter equation. 
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The solution is given by 



K\k + 1 



u 



K [k-l 



[c + u\[k + c — u\ 
[c — u][k + c + u] 



(2.11) 
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Here the constant c is arbitrary. We choose to restrict c to he in one of two separate regions 



parameterised by A) x 
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-1 < 6 < 1; B) X 



2c __ ^2b 



satisfy < u < |6| < 1 in both regions. The overall factor in i^' it 1 



Graphically, we represent Kla u]hy 



T < 6 < 1. We further restrict u to 
k 



k 



u is determined later. 




2.4 The boundary transfer matrix. 

We define the boundary transfer matrix 



The boundary Yang-Baxter equation implies [1^, |S|] 

[T;}\u),TPiv)]=0. 
Graphically, (u) is represented by the following half-infinite transfer matrix 
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(k) 

We also define the boundary Hamiltonian ' in terms of the transfer matrix as 



H 



2x 



dC 



lC=i- 



(2.14) 



Our first question is about the ground state configuration, i.e., the largest eigenvalue eigenvector 
of T^^^ (u) in the limit x ^ 0. We know that in the bulk theory the ground state configuration 
of heights is given by the alternating sequence of / and I + 1. Because we have chosen to fix the 
heights at the right-most corners to be equal to k, it follows that the possible choices are I = k oi 

( k \ 

k — 1. Let us compare K IkzLl ^ u j. In fact, in the limit x — > we have 



K{k + 1 



K\k-l 
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> 1 if 6 > 0; 
< 1 if 6 < 0. 



(2.15) 



This is true in both regions A and B. Therefore, the ground state configuration for 6 > is given 
hj I = k^ and for 6 < it is given by I = fc — 1. 

Next, let us fix the overall factor in the boundary weights. The normalisation will differ de- 
pending on whether 6^0, and we shall use the corresponding notation iC^^ {ji+e^ uj for the 

normalised boundary weights. 



First, suppose 6 > 0, and set 
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The corresponding boundary transfer matrix T^^ {u) satisfies 



5 y-u)f'i\u) = l, 

(fc) _ [2-2u][c + u][k + c- 



R{2u-\)[2u][l + c-u][k + c-\ + 



(2.16) 
(2.17) 



Since {T^\u)^ is a commuting family, the eigenvalues of T'^'{u) also satisfy these equations. As 



suming the analyticity of logi^^\u) for the largest eigenvalue i^^\u) in a neighbourhood of the an 
nulus X < Id < 1, we determine t^^\u). After that we choose the normahsation of K^^ 

in such a way that Tj^\u) has the largest eigenvalue 1. The result is as follows 
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/(C) 

This normalisation of K 
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uj is valid for both regions A and when 6 > 0. The case 
6 < is similar. It is even not necessary to repeat the calculation because we have a symmetry 
(c, k, lb) — > (— c, r — /c, =p). Therefore, we have 
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In conclusion, we have fixed the normalisation of the boundary weights for 6 > and 6 < 
separately so that the largest eigenvalue of the boundary transfer matrix is 1. 

For completeness, we give two relations for T^'^ (u) and the related relations for K^'^l k+e ^ u ] 
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3 Bosonization 

We follow to obtain a bosonic realization of the spaces Ci k, their duals, and the vertex operators. 



3.1 Boson Fock spaces 

We use the bosonic oscillators 



[Pm ) Pn] 

[P,Q] 



f^^ — n i Om+n,0, 



(3.26) 
(3.27) 



where [k]x = ■ The Fock spaces are 



where 



Fi,k = C[/3_i,/3_2,---]|^fc), (3.28) 
J^^j, = (/,A;|C[A,/?2,---], (3.29) 
{Uk\Uk) = 1, (3.30) 



I3m\hk) = Oifm>l, (3.31) 
{Uk\Pm = Oifm<-l. (3.32) 



The operator P acts as 



r — 1 



y 2(r - 1) ~ ^ V ^ "^''^^ ^"""^ ^^"^^^ 

We also use operators K and L which act on JF; fc, JF^*^, as k and Z, respectively. We use the 
convention of left action on J^i^k and right action on J'^f^- Namely, the composition AB of operators 
A and B is such that B acts first on J^i^k and A acts first on J^ij^. 

3.2 Type I and II vertex operators 

We define the type I vertex operators 



r-l 



'r^(Q+iPlogO-Y: 



$+(C) = C~ ■■ e :, (3.34) 

- -'i^*^<C)-^(^)^^^[^^. (3.-) 



2(r-l) 



A(^) = z— : e -"^o :, (3.36) 



where 



( = x^", z = x^", (3.37) 

\ _ ^^Ct^'— 2 . ^2r-\ Joo^.^ ,X ,X Joo 

A-X r-l^X ,X Jcx)/ 2.^,4 ^2r^ /'™2+2r. ^,4 ^2r^ " 

The contour C for the z-intcgration in $-(C) is chosen in such a way that the poles from the factor 
[u + V — ^] at z = x^'^'^"'^ /C are inside if n G Z>o (outside if n € Z<o))) and the poles from the 
normal ordering of the product $+(C)^(z) at z = l/x^~^^"'^( (n G Z>o) are outside. The operator 
^±(0 acts as 

^±(0 : :Fi,k ^ J^i,k±i, (3.39) 
^±(0 : n,k±i^n,k- (3.40) 

The value of K in $-(C) is k in both cases. 
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The type II vertex operators are 



^-(0 = i.TT— n(CM^) "r 1]/ ^ (3-42) 

(3.43) 



where 



«m = (-ir ., ^ /3m, (3.44) 

[(r - l)m]a 



2r 



^(Q-iPlogz)- ^ ^(x^+x-™)^-™ 

B{z) = -.e"' " :, (3.45) 

[u]' = x^-"e,2(.-i)(2;2"). (3.46) 

The contour C' for the z- integration in *I'1(C) is chosen in such a way that the poles from the 
factor [ti + V + at z = are inside if n G Z>o (outside if n G Z<o) and the poles from 

the normal ordering of the product ^*j^{C,)B{z) at z = (n G Z>o) are outside. The 

operator ^'±(C) acts as 

^liC) : J'i,k^J'l±i,k, (3.47) 
^i(C) : ^r±i,fc--^z:fc- (3.48) 

The value of L in ^'l(C) is I in both cases. 

The vertex operators satisfy the following commutation relations. 

<f..(C2)f.,(Cl) 

E^:(C)^e(C) = l, (3.50) 

e 

$,,(C)C(C) = 5e'e, (3.51) 

where 

m) = ^f-eix\). (3.52) 

Comparing (p9|)-(^) with (l2?7|)-(|2^, we see that ^eiC) and <I>*(C ^) provide a Fock space 
realisation of the half-infinite transfer matrices $('^+<^''^)(C) and {'=.*=+'^)(^-i) respectively. 
We also have 
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u I is the Boltzmann weight with r replaced by r — 1, and 
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Our bosonization is slightly different from than that of reference Q. Firstly, we have changed 
C to in the definition of ^>±. Secondly, our involves an extra factor of (—1)^"^ over the 
corresponding object ^~ of 0. This factor arises because our Boltzmann weights W differ from U 
of by a — sign in (^). This difference is simply a gauge transformation. Finally, our am is equal 
to (— l)™'am of 0. We introduce this factor in order to give the correct sign in the commutation 
relations ( |3.54| ). 

3.3 Screening operators 

We define two screening operators X and Y that have the properties |T^ . 



They have the following bosonization 

X = 
Y = 



y ■ ^i,k ^i,k-2 {^ik~2 ^ik)^ 

[ci>,(c),x] = o, [^:(c),r] = o, [x,y] = o, 

0, Y'' = 0. 



dz 
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It is proved in that these expressions obey ( |3.62| ). 

Fix / (1 < / < r — 1) and k (1 < / < r). We have the following complexes. 

X'2=X^ -rr X_i=X'-i-' Xo=X' Xi=X'-i-' 

••• ^ ■^2{r-l)-l,k •^l,k ^ •^-l,k " ■■■ [6.bb) 

••• -'2{r-l)-l,k " ■^l,k " ••• (6.0 t) 



^l,2r-k — " ^l,k ^ ^l,-k ^ • • • (3.68) 



■ ■ ■ /,2r-fc l,k " -> l-k " ■ ■ ■ 

We assume the following cohomological properties. 

KeiXj/lmXj^i = Oifj^O, (3.70) 

KerX*_i/ImX* = if j / 0, (3.71) 

KerXo = Ker^o, (3.72) 

ImX_i = Imy_i, (3.73) 

KerYj/lmYj^i = if j / 0, (3.74) 

Kery/_i/Imy/ = if j / 0, (3.75) 

KerX*i = Kery!i, (3.76) 



We will make the identification, 



and conjecture that the coupling 
induced from the coupling 



ImXo* = Imyo*- (3.77) 

= KerXo/ImX_i, (3.78) 

A% = KerXli/ImX*, (3.79) 

CI, X ^ C, (3.80) 

Tr, X - C (3.81) 



is non-degenerate. 

In summary, we have bosonized the half-infinite transfer matrices (I)('^^^''^)((^) and $*('^''^^i)((^~i) 
by means of the type I vertex operators <5±(C) and $^(C~^), and also introduced type II vertex 
operators, which will play the role of creation operators of particles as we will see in the following 
section. 
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4 The Boundary States 



4.1 Construction in the Fock space 

In Section 2, we chose the normahsation of the bulk and boundary weights such that the maximal 
eigenvalues of the bulk and boundary transfer matrices in the large lattice limit were 1. In this 
section we describe how to construct the corresponding eigenvector in the bosonic Fock space T^.k 
or J^fc-i,fc) depending on whether 6 > or 6 < 0. We will then show that the eigenvector so 
constructed actually determines a non-zero residue class in the corresponding subquotient Ck^k or 
i^k-i,k- 

The maximal eigenvectors fc)^ and |A; — 1, A;)^ are defined in both regions A and B (defined 
in Section |2!^) by, 

^^^(Ol^,^:^^ = \k,k)% e J'k^k, for6>0, (4.1) 
T'i\0\k - l,k)% = \k-l,k)%eJ'k-i,k, for6<0. (4.2) 
We make the Ansatz that 

\k-i,k)% = e^^'''\k-i,k)M i = 0,l, (4.3) 

where 

, (4.4) 

2 '—^ m\rn\x\{r — l)m\x m 

m>0 ^ J-^LV J 1-1. m>0 

The adjoint action of e ' on the bosonic oscillator modes is that of a Bogoliubov transformation. 
Namely, 

p-Ff'' R pK''' - R -R I [mh\{r-l)m]^ jjc,k 
e p^e « - Prn P-m^ [2m], [rm], ^m,i^ 

where m > 0. The coefficients D'^^ are determined by solving conditions ( [4.1| ) and ([4.2|). Rewriting 
these conditions in terms of the bosonic vertex operator expression for T^\c,) ( |2.12| ), and using 
the inversion property ( |3.51| ), we obtain 

k 



K^'^ { k+e 



k 



n $e(C)e^o \k,k) = <^e{C )e''o \k,k) for 6 > 0, (4.6) 



k. 

F^-^iu 1 ;,\ _ ^ /a-1\„-F; 



n U>e(C)e^i \k-l,k) = (^.{Cle^' \k-l,k) for 6 < 0. (4.7) 



Solving these conditions for e = +, we find 



r.c,k _ [{k-l)mU{r-2c-k)m]+ ^ ( [m/2Urm/2]i \ 

~ [{r-l)ml '™ 1^ [(r - l)m/2] J ' ^^"^^ 

k _ [r + l-k)mU{2c + k)m]t ^ f [m/2Urm/2]t \ . . 

~ [{r-l)ml 1^ [(r - l)m/2] J ^^■^> 

j^c,k , [rm]x[{2c + l)m] + 

~ [(r-l)m]. ' ^^-'^^ 
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where 

[/c]+ = x'^ + x^^, 

X if m is even; (^-H) 
otherwise. 

In order to show that ( |4.6| ) and ( |4.7D hold for e = — , with D"^- given by ( |4.8D and ( [4.9D , we first 
consider (|4.6| ) for = 1. We find that, 

^y2-Kizy^ ^ ) (xC^;x^'-)ooK/2;x^'-)oo ^ |J-,-l;-u- 

(4.12) 

The integral vanishes because of the anti-symmetry of the integrand under the change of variables 
z z^^. The proof that ( |4.6D and ( [4.7| ) are valid for e = — and general k is inductive and is given 
in Appendix |^. 

In a similar manner we construct dual maximal eigenstates in the Fock spaces J^^ and 
These states are determined by the requirements, 

%{k,k\T'i\(:) = %{k,k\eH,k for6>0, (4.13) 

%{k-l,k\Ti,'\c) = %{k-l,k\en-i,k for6<0. (4.14) 
As above we make the Ansatz that 

%{k-i,k\ = {k-i,k\e^'^'\ (4.15) 

where 

^c,k ^ _1 „4m [2m]x[rm]a, ^2 , "^m,* 
9 

m>0 ""[""jaJLV -^/""ja; m>0 

e i now produces the Bogoliubov transformation 

^_G^''= _ « ^4m« , N[(r - l) m] ^c,fc 

[2m] [rm 



m[m]x[(r- l)m]a; f-'^ m 



= /?-^-x^-/3^+ ^";^^^:,,_;'-^ i?:;^, (4.17) 



e^' /3me-^« = (3m, (4.18) 

where m > 0. In terms of bosonized vertex operators, ( [4. 131) and ( |4.14| ) become 

k 



%{kMm-')K^>^[k+e^ 



%{k-lMm-')K'iHk+e^^ 



= |(A:,fc|cI>:(C) for fe>0, (4.19) 
%{k-l,k\^l{C) for 6<0. (4.20) 



Solving these equation, we find 



^-•0 - ""^ [(r - l)m] I [(r - l)m/2], ' ' ^^'^^^ 



^c,k _ 2m [ir-k-l)m][{2c + k)m]t omn ( [m/2Urm/2]t \ 

[(r - l)m\ \ [{r - l)m/2\^ I 



:^c,k , ^2m [r m]^[{'^ - 2c)m ] + 
[(r — l)m] 



^;^fe^+^2n. i--mv- (4.23) 
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4.2 Boundary vacuum states in the subquotients 

We have constructed maximal eigenvectors \k — i, k, A;)^ € !Fk-i,k, and dual eigenvectors ^^{k — i, k, k\ S 

^. In order to establish that these vectors actually give rise to the eigenvectors in the subquo- 
tients Ck-i^k = KerXo/ImX_i and C^^^j^ = Ker X'^i/IiuXq respectively, it is sufficient to show 
that, 



\k-i,k)% G KerXo = Keryo C^fc-,,fc, (4.24) 



"".{k-^kl G KerX*_^=KerYl^cn-^,k, (4.25) 
''^{k-i,k\k-i,k)% / 0. (4.26) 



If these three conditions are true then it is simple to show that 



B 



{k - i, k\P\k - i, k)% = %{k-i, k\P\k - i, k)%, (4.27) 



where \k — i,k)'^Q is the residue class of \k — i,k)'^^ in Ck-i,k, %{k — i,k\ is the residue class of 
^(A: — z, k\ in C*^_i and P is an operator that commutes with the coboundary operators X and 
Y . Local operators of the of the theory are expressed in terms of products of vertex operators 
<I>e(z), and thus an example of an operator P. Hence if ( 4.24| ) - ( 4.25| ) hold we can calculate local 



correlation functions of the theory as matrix elements of the form ^(/c — i, k\P\k — i, /c)^. Form 
factors are also in this class because they are the matrix elements of products of type I and type II 
vertex operators which commute with X and Y respectively - see ( p. 61 ). 



We establish condition ( [4.24 ) as follows. For |1,1)^ = e^o' |1,1), one can show the result 
directly. Namely, 

"-m [m] , 

e-^o' Xe^o' |l,l) oc i — — (z-z-i)e™>o |-l,l)=o. (4.28) 

J 27riz 

This integral is zero because of the anti-symmetry of the integrand under z z^^. For general 
I A; — i,k)'^ this technique is rather complicated and it is simpler to proceed by other means. We 
observe that different boundary states are related to each other by the action of vertex operators. 
In particular we find that 

^+{x^')\k,k)% ~ \k,k + l)%-\ (4.29) 
Res^^_^i+2c^X{(,)\k-l,k)% ~ \k,k)%. (4.30) 



Here ~ means equal up to scalar functions. We shall discuss ( [4.30 ) and other relations in more 



detail in Section 6. These relations are useful in this context because the screening charge X 
commutes with ^{(,), and Y commutes with 'I'*(^). Thus 

X|l,2)^ ~ cI>+(x2(^+i))X|l, 1)^+1 = 0, 

Y^\2,2)% ~ Res^=_,i+2c^;(Ol^'|l,2)^ =0, 

where |1, 1)^^ is an analytic continuation in c of |1, 1)^. By repeatedly applying vertex operators to 
construct the general state \k — i, k)'^, we can show that \k — i, k)'^ S KerXo = Ker Iq as required. 
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We establish condition ( 4.25| ) in a completely analogous manner. First we show the result for 
^(r — 2, r — 1| directly using, 



(r — 2, r — l|e*^i' Xe oc (r, r — 1| 



—— (zx -z X )e 

I-KIZ 



(4.32) 

1^-2 



This integral vanishes because of anti-symmetry under the changes of variables zx^ ^ z "x 
Again, different dual states are related to each other under the action of vertex operators. We find. 



B 



{k-l,k\<^*_{x'^) ~ ^ 



Res^=_,i-2c^(fe,A:|M/_(r') ~ ^(fc-l,A:|. 



(4.33) 



(4.34) 



Thus 

%{r-2,r-2\X ~ ^"^(r - 2, r - l|X$l(x2(^-i)) = 
^(r-3,r-2|y2 ^ Res^=_3,i-2c ^(r - 2, r - 2|y2$_(^-i) = etc. 

The general result '^{k — i,k\ G KerX^^^ = Ker y!';^ follows. 

Finally, we show ( 4.26| ) by computing ^(A; — i,k\k — i,k)'^^ explicitly. We do this by using a 
decomposition of the identity on J^i k which employs coherent states. Define the coherent states 



^ g„oo (fl = (/, A;| e'">o , where. 



1 [2m]j,[rm.]3: 
m [m]x[{r—l)m]x ' 



(4.35) 



It is simple to show that 



(4.36) 



Its is possible to decompose the identity on J^i ^ as 



This is shown using the Gaussian integral 



m>0 



(4.37) 



/ n d^rnd^m^m(^ rnX) m>0 
— oo m>0 

= n (-det(^„))-2e ™>o 

m>0 

Inserting this decomposition of the identity between e ^ e » and using (4.38), we find. 



(4.38) 



{k — i,k\k — i,k)'^ = [k — i,k\e^ e » \k — i,k) 



1 



1 



^exp I J2 



^4™(^c,fc.2_22:,<:.fc£;=.fc+(£;=>'=.)2 



(4.39) 
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For i = this is equal to 



(a;2+4r-4c-2fc .a,4^^2r)^(a,2+4c+2fe . a,4^a,2r )^ (a,2+2r .^4^^2r)2^ (a,2r+2-2fe . a,4^^2(r-l))^ (^2r+2+2fe . 2.4^a,2{r-l))^ ^4.4Uj 
/^4r+4-4k-4c . ^8.^2r\ (-„4+4c . „8 .™2r -i /'„4r+4-4c . „8.„2r-i /'„4+4fe+4c . „8 . „2r -i 

^ ■p7+2=4l^rfcT^X;;2(?^rTyj^^^ 

From the exphcit form of D'^^ and E^-, we see we can recover the i = 1 normahsation from the 
1 = result by the transformation (c, A;) {c + r,k — r). 



5 General Eigenstates and The Boundary S-Matrix 

5.1 General Eigenstates of T^''^(C) 

We construct general eigenstates of Tg^-* {() by using type II vertex operators. The technique closely 
parallels that used for the XXZ model |, |]. From (^J^) and (|^) it follows that 

Tif\OmO = t{C,OmOTl}\C), where, 

(5.1) 

t{C,0 = x(C/C)x(eC/^') = x(C/C)x(i/(CC)). 

From the definition ( |2.14| ) of the Hamiltonian and the property r^'^^(l) = id (which follows 
from the definition ( p.l2| )), it also follows that, 

= e(0^:(e), where, 

(5.2) 

6(0 = f sinh(^)dn(2I^), e = -e2^^ 

Here, dn is the Jacobi elliptic function with nome x, and / and /' are the associated half-period 
magnitudes [Q]. 

Thus it is possible to construct general continuous eigenstates of the form 

nA^i)---nA^n)\k,k)% e J^k+e,+-+s„,k for6>0, 

^l,(.(l)---^lSn)\k-l,k)% G J^k-l+e,+-+er.,k for 6 < 0, , 

with = 1. The physical picture of these excitations is of n kinks or domain walls that are free, 
i.e., not bound to the right boundary. 



5.2 The Boundary S-Matrix 

The boundary S-matrix describes the scattering of these kink states off the right boundary. Fol- 
lowing [|lC| ], we define the boundary S-matrix mJ'^ *'^^(^;r, c) by 

^m\k-i,k)% = M('-''''\a;r,c)^t{r')\k-i,k)% fori = 0,l, (5.4) 

%{k-i,k\^,ir') = %{k-i,k\^,iOM^^-''''\a;r,c), (5.5) 
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a], (5.6) 



where ^ = Solving condition and ( |5.5D We find, 

where 5 = jj^, and K' indicates that r is replaced by r — 1 in the boundary weights. The 
boundary S-matrix elements are thus given in terms of the boundary Boltzmann weights ( |2.18| )- 
( |2.21[ ) analytically continued to shifted arguments. Again, this is very similar to the situation for 
the boundary XXZ model The sole effect of the 6 is to produce the minus sign in the shift 

In order to discover discrete eigenstates of T^\(^) other than the maximal ones, we should 
examine the boundary S-matrix elements in search of simple poles in the physical strip. We 
identify the boundary physical strip as x < |^| < 1, and find just two simple boundary S-matrix 
poles which can lie in this region for the allowed values of c. Firstly, we find that there is a pole in 
M^'''\a; r, c) at ^ = —x^~'^'^. By explicit calculation, we have shown that this arises from a simple 
pole ~ l/(^ + x^"^"^) in ^*_{^)\k, k)% with residue 

ReSi:=_,i-2c^*_iC)\k,k)% ~ \k-l,k)%,. (5.8) 



where \k — l,k)'^ is an analytic continuation of ( [4.3| ) in c. Equivalently, this boundary S-matrix 
pole arises from a pole in %{k, A;|^'_(^~^) at ^ = —x^~'^'^. The residue at this pole is 

Res^=_,i-2c%{k,k\^-{r')-'B{k-l,k\. (5.9) 

The second pole that can lie in the boundary physical strip occurs in 
_^i+2c_ pqIq arise from one in ^^(^)|A; — 1, /c)^ with the following residue, 

Res5=_,i+2c^;(0|fc-l,A;)^ ~ \k,k)%,, (5.10) 

or equivalently from 

Resg=_,i+2c %{k - 1, A:|^+(r') ~ iik, k\. (5.11) 

There are no poles in M^'''\a;r, c) or ^'''\a;r, c) in the physical strip for either of the 
parameter regions A or B. 

We shall now discuss the regions A and B separately and in more detail. First consider regime 
A, in which x^^ = — x^^, —1 < 6 < 1. For 6 > there is a unique ground state with energy 
eigenvalue zero. This is |A;,A;)^. For < b < ^ the pole (|5.8| ) in M^'''\a;r,c) at ^ = x^^^^ lies 
in the physical strip. The state \k — 1, k)'^ can be interpreted as a bound state of a '— kink' with 
the boundary. The energy of this bound state increases from zero at 6 = up to a maximum value 
of e(l) at b = ^. This maximum value is equal to the minimum of the energy band of a single 
free kink state (i.e. one of the form ^*_{^)\k, k)'^^ with |^| = 1). For b > ^ the pole moves off the 
physical strip, and the state \k — 1, A;)^ can decay by emitting a kink. For 6 < the unique ground 
state with energy zero is \k — 1, k)'^. The pole in ^'''\a; r, c) at ^ = x^'^'^'^ lies in the physical 
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strip for — ^ < 6 < 0. In this region \k, A;)^ is a bound state of \k — 1, k)'^ with a + kink. Again, 
for b < the pole moves off the physical strip and |A;, fc)^ is unstable against emission of a free 
kink. 

Now consider regime B, which is parameterised by x^'^ = x^^, — 1 < 6 < 1. For < 6 < ^ 
the pole (5.8) in M^'^\a;r,c) at ^ = —x^~'^^ still lies in the physical strip, and the energy of the 
bound state | A; — 1, A;) ^ is outside the energy band of a free kink. However, in this case, the bound 
state energy is greater than the maximum energy e(— 1) of a free kink. As b increases the energy 
of the bound state \k, A;)^ decreases and reaches e(— 1) when b is equal to the critical value |. For 
6 > ^ the pole moves off the physical strip and \k — 1, A;)^ can decay by emission of a kink. This 
mechanism for the stability of a bound state in a lattice model was also observed for the XXZ 
model in [^]. For 6 < 0, again, and for similar reasons, there is a stable bound state |A;,A:)^ when 
-i < 6 < 0. 

We remark that the normalisation ( [4.39|) also confirms the above discussion, because %{k, k\k, k)"^ 
has has a pole at x'^'^ = ±x~^, and ^(A: — 1, A;| A: — 1, k)^ has a pole at x'^'^ = itx. 



6 The Scaling Limit 

In this section we consider the x ^ 1 scaling limit of the bulk and boundary S-matrices. Param- 



eterising ^ = x^", we shall identify nia with the usual rapidity parameter 9 in this limit. Our 
boundary physical strip becomes < Im^ < 7r/2 as in The pole at = x^~^l^l in region A 

thus lies on the imaginary axis in the 9 plane. It is a physical strip pole that may be interpreted as 



a bound state. The pole at ^ 



„l-2|f)| 



in region B moves off to infinity in the scaling limit. The 



energy of the associated state becomes infinite, as the one particle energy band becomes infinity 
high, and it should no longer be considered as a bound state. It is for region A that we shall 
construct the scaling limit of the bulk and boundary S-matrices. 



6.1 A gauge transformation 

Before taking the scaling limit we carry out the following gauge transformation: 

/ [7714]' [7712]' I [mi] 



S 



?7i4 mi 

7713 7712 



a 



V [7771]' [7773]' 



[7773]' 



w 



7774 7771 
7773 7772 



a 



(6.1) 



a 



[k-i + e]' 
[k - i]' 



(6.2) 



[l + e]' 
[I]' 



a/2 



\[l + e]' 



(6.3) 



(6.4) 



^In the boundary physical strip was at one point incorrectly identified as < Im^ < tt. This was corrected 
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where ^ = x^"^. The effect of this gauge transformation is to yield bulk and boundary S-matrices 
S and M that have the conventional crossing and boundary crossing/unit arity relations (a rather 
similar gauge transformation was carried out in reference in order to yield a crossing symmetric 
S-matrix from RSOS Boltzmann weights). Written in terms of these gauge transformed objects, 
equations ( p.53 ),(5.4) and ( ^.SD become, 



*:,(e2)^:,(6) 



e[+e'2=ei+e2 



L + ei+S2 L + e'2 
L + ei L 



Q2-ai (6.5) 



%{k-i,k\^e{r' 



M^''-''^\a;r,c)^*{C^)\k-i,k)% for i = 0,l,, 
%{k-i,k\^,{OM^'^-''''Ha;r,c). 



(6.6) 
(6.7) 



The S-matrix (6.1) obeys ( |2.4| ) and ( |2.5D with the change W ^ S. The crossing relation is modified 
to the more conventional 

a d 



S 





c 




r 




a 


V a 


b 





s 



b c 



1 - a . 



(6.^ 



The properties of the boundary S-matrix M can either be read off from the explicit expressions 
( ^I^ ) and (|5.7|) , or derived using the definition ( |6.6D and properties (|6.5| ) and (p. 57 ) of vertex 
operators. They are: 



£'{+e'2=e\+e2 



s'(+e'2=ei+£2 \ ' + £2 



l + ei + £2 1 + e'{ 
1 + 62 I 

l + ei + £2 1 + e'i 



I 



ai — 02 j S 
ai + 02] S 



1 + 81+62 1 + £2 

I + e'i I 

I + £1+62 l + e'2 

I + e'i I 



ai + a2 mS';'^(6)mS'''^(6) 



a,-a2]Mg'\i2)Mf,;'\i,), 



where ei, e'l and e'2 are fixed with ei + £2 = e'l + ^21 and I = k or k — 1, 

0'''\c)MpHr^) = 1, 



e' \ l + e I 



(6.9) 



(6.10) 



6.2 Scaling behaviour 

The S-matrices are constructed in terms of functions that have the following behaviour in the limit 
p — > 1: 



r(z) 

Gp(p^) sin(7r2;) 
0p(P^') sin(7rz') 



(6.11) 
(6.12) 



We shall consider the two cases r = 4 and r = 5, separately, and in detail. The scaling limit of the 
r = 4 model should correspond to the massive <A(i,3) perturbation of the critical Ising model [ [To| . 
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That of r = 5 should correspond to that of the massive (p(i^3) perturbation of the tri-critical Ising 
model m, 0]. 



6.3 The r = 4 model 



s 



s 



2 1 

1 2 

1 2 

2 1 



02 - ai ^*+{(iWM2), on ^2,fc 



02 - ai on J^i^k, 



where 



S 



2 1 
1 2 



a2 — oi 



5 



1 2 

2 1 



02 — «! 



-1. 



For r = 4, there are just two S-matrix elements. Reading off from equations (|6.5| ), (6.1), and 
(0)-(|2l), we find 



(6.13) 



(6.14) 



(6.15) 



This is the behaviour expected for a free fermion theory. In order to define the scaling limit 



of the boundary S-matrix elements given by (|6.2| ), (|5.6| ) and (^j), we must scale ^ 
consider region A in which we have the parameterisation x^'^ = 
that in the limit x — > 1 



X 



2a 



We 



x^^ with -1 < 6 < 1. We find 



M^''^\a; r = 4, c) = M^L'-^^a; r = 4, -c) 



(2,3), 



tan(7r(a/2 - 1/4)),, 



(6.16) 



M 



(2,2) 



(a; r = 4, c) = M|^'^^(a; r = 4, -c) 



, , , ,,,, /cos(7r6) + sin(7ra) , , 
tan(7r(a/2 - 1/4)) —- — ( (6.17) 



,cos(7r6) — sin(7ra) 

We now compare these expressions with the boundary S-matrices for the Ising model derived in 
reference [pX|]. We find that ( |6.16D and ( |6.17] ) agree with the boundary S-matrices for 'fixed' and 
'boundary magnetic field' boundary conditions given in (4.10) and (4.27) of |1C| if we identify the 
rapidity 9 = iria, and the magnetic field parameter k = cos(7r6). In addition, setting 6 = in ( |6.17 ) 



gives the boundary S-matrix for 'free' boundary conditions as given in (4.16) of [10|. 
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6.4 The r = 5 model 



There are just four independent non-zero S-matrix elements in this case. Their scahng behaviour 
is given by, 



i 

where, 
F(a] 



3 2 

2 1 

2 1 

1 2 

2 1 

3 2 



a 



a 



a 



a 



S 



S 



S 



S 



3 2 

2 3 

1 2 

2 3 

2 3 

3 2 

2 3 

1 2 



^2^ "cos(f (a- l))F(a) 



a 



a 



a 



a 



V2^ "cos(f(a + l))F(a) 



-^/2"cos(fa)F( 



-^/2" sin(f a)F(a) 



r(n+i)r(n 



)r(n + f 



n>or(n-i)r(n + f)r(n-f + i) 



.ttq;, 
cos(— ) 



1 



exp 



-Boi9), 



-Biie), 



dt sinh(at) 



(6.18) 



4 Jo t cosh^(t/2) 



(6.19) 



The function F{a) has the properties 

F{a)F{-a) = l and F(l - a) = F(a). (6.20) 
Bq{6), Bi{6), ^o(^)) ^i(^) the S-matrix elements given in references |16, (note however 



that F{a) is convergent, unlike S{d) of references [16, |TT[). The representation of F{a) in terms 
of r functions was obtained using the limiting behaviour of {p^',p)oo discussed in Section 7.2. 
The integral representation was constructed by taking the limit of the infinite sum that occurs in 



\og{p^;p), and is equal to the expression given in reference |17|. It is simple to show that the two 
representations agree by making use of the integral representation of log(r). 



The symmetry of mJ*^ *'*^''(^; r, c) 



1 k+i,r means that there are four inde- 



pendent boundary S-matrices. They have the following scaling limits: 

^"tan^d -i)P^in(a), 



M'^'^\a;r = 5,c) 



Mi^'^)(a;r = 5, -c) 



M 



(3,3) 



M. 



(2,2) 



(a; r = 5, c) 



(a; r = 5, c) 



M^^'^\a;r = 5,-c) 



e^^tanh(f - f)i?(_i)(0), 

^/9"tnn7r('° ^ P ■ ( ( cos(7rb)+sin(7ra) 

e^^aanh(f - |)P-(^), 



(2,3) 



(a; r = 5, — c) 



V2-"tanvr(f-l)P_(a)(H2|g 



fe)+sin(7ro) 
fe)— sin(7ro) 



F{a-b 



-i)F(a + 6+i)(^cos(vr(| + i)Tsin(^; 
eT^aanh(f - ^)R±{9). 



(6.21) 
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where 



„ r2(n-f)r(n + f + l)r(n + f -i) _ /_ foo dt smh(at) smh\t/A) ^ 



(6.22) 



„>o r2(n + f )r(n - f + l)r(n - f - i) "V -^0 * sinh^ W 
-Pmm(a) has the properties 

-Pmm(a)-Pmm(-a) = 1, and Pmin{l - a) = \/2sin(7rQ/2)F(2a - l)Pmin{a). (6.23) 

The latter is shown most simply by using the integral representation of both sides. R{-i), P-{Q) 
and R±{6) are the boundary S-matrices of reference [0], and 7 is defined in ll| as e^'^*'^ = 2 (we 



identify a = and irb = ^ — ^ , where ^ is the magnetic field variable of reference |Tj 

We have checked that the scaling limits of our bulk and boundary S-matrices obey the continuum 
Yang-Baxter, crossing, and unitarity conditions, and we should comment on overall factors by 
which our results differ from those of references |16, |l^]. If we wished, we could remove the overall 



— sign in our S-matrices by a gauge transformation. However we choose to keep it because it 
introduces the tan7r(a/2 — 1/4) factor in the boundary S-matrices. As pointed out in reference 
|11], we expect a simple pole in M^'^\a;r,c) = M^''^\a;r,—c) at q = 1/2 when 6 = 0. The 



'CDD factor' (cos(7r6) -|- sin(7ra))/(cos(7r6) -|- sin(7ra)) however has a double pole at a = 1/2 when 
6 = 0. So we need the factor tan7r(a/2 — 1/4) in order to turn the double pole into a simple 
pole. This situation is just the same as for the r = 4 Ising model where the same factors are 
present in M^''^\a; r,c) = M +(-^'^) (a; r, — c). The absence of the remaining e'^^ in the boundary 
S-matrix elements of [|ll| appears to be a simple error. Its presence is necessary to ensure boundary 
crossing/ unitarity. 



7 Discussion 

In summary, we have extended the bosonization scheme of reference to the boundary ABF 
models. Specifically, by making use of bozonized vertex operators we have constructed boundary 
states \k — i, k)^, and the associated boundary S-matrices. We have initially constructed \k — i, k)'^ 
in the bosonic Fock space J-k-i,ki and then shown that it determines a non-zero residue class in 
Ck-i,k under some assumptions on the cohomological construction proposed by Lukyanov and 
Pugai. For 6 > 0, \k,k)'^ is the vacuum state, and \k — l,k)'^ ~ Res^=_x.i-2c ^^1(^)1^;, A;)^ is 
interpreted as a bound state of \k, k)'^ with a — kink. This bound state is stable for < 6 < i. 
Conversely, for 6 < 0, |A; — 1, fc)^ is the vacuum state, and \k, k)% ~ Res ^=_a^.i+2c ^'+(^)|A; — 1, fc)^ 
may be interpreted as a bound state of |fc — 1, fe)^ with a -|- kink. This bound state is stable for 
— ^ < 6 < 0. General continous eigenstates of the boundary transfer matrix are constructed by 
acting with the type II vertex operators ^* on the discrete states \k — i,k)'^. 

The scaling limit of our bulk and boundary S-matrices should be those associated with the 
^(1,3) perturbation of the c = 1 — 6/r(r — 1) rational conformal field theories. For r = 4 and r = 5 
we have checked our scaling limit against the results for these perturbed boundary conformal field 
theories which are available in the literature. 
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There are several directions in which this work may be exploited and extended. Using our 
bosonized expressions for the bounds states and vertex operators, it is possible to write down 
integral expression for arbitrary equal time correlation functions of local operators of the theory. 
This may be done by a simple extension of the technique presented in Q. Integral expressions for 
form factors may be derived in the same manner. Difference equations for the correlation functions 
and form factors may be derived following the approach of The results could be extended to 
^n-i ^^^^ models using the bosonization scheme of reference [19|. 
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A Normal ordering relations and summation identities 
A.l Normal ordering relations for vertex operators 



We list the normal ordering relations for the vertex operators discussed in Section 3. We abbreviate 
AB =:AB:C to AB =: : C, and use C = a;^" and z 



$+(Cl)$+(C2 

A{zi)A{z2 
A{zi)A{z2 

^;(Ci)*;(C2 

B{z)^XiC 

$+(Ci)^';(C2 
*;(Ci)^+(C2 

^>+(Ci)^;(C2 

^+iC)Biz 

n(OA{z 
nmiz 

A{z,)B{C2 
B{z2)A{Ci 

A{zi)B{C2 

B{zi)B{z2 
B{zi)B{z2 



■ ■ ^1 {x''-QlK2\x'',X^^)oo{xHl/C.2;x'',X^^)oo ■ 

■ ■ {zCx;x-^noo ' 
. . (Z^) {{zO-^x^-^;x^n^ 



A(^)$+(C) 



l^+u+v] ' 



■.zr^il-Z2/z^)j-&^f^^ 



[■i;i-t'2— 1] 

^1-1^2 



^Aiz2)Aiz,), 



. >-(2(;£^) (Cl/C2;x^a:^('-l))oo(x2('■+l)Cl/C2;x^x^('-l))o 
■ ^1 (x2'-Ci/C2;x4,x2('-i))„c(a;2Ci/C2;a;4,a;2(r-i))<^ 



■■ c 



(2Cx-l;a;2{''-i))oo ' 
{izO-^x-^;x^n 



n(c)^(^)}pSi7, 

. . A {-xHi/C2;x^)^ 
■ ■ ''I (-a;Ci/C2;x4)«, ' 

. . /-^ (-x3Cl/C2;X^)cx. 

• • ''I (-<i/C2;2^'*)oo ' 

(Ci/C2)^?^^{^n(C2)<&+(Ci), 

: :z(l + (Cz)-i), 
: :z(l + (C^)-i), 

^(^)n(c), 



2^(l+a;22/2i)(l+a;-i22/2i) ' 

. . 1 

z'j{l+xzi/z2){l+x-''-zi/z2l) '■ 

B{Z2)A{C,), 

2r 

^^^B{z2)B{z,). 



■ ' (1 - ^'il^^)%2rll'll];{'^2(r-l)]l 



(A.l) 
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A. 2 Summation identities 



E 

gm>0 
gm>0 



m [27n]x[(r—l)m 



[m]j,[(r-l)m: 



[m] 



E 

m [rm]a; 



gm>0 



E 



m [(r— l)m 



(2;x4,x2('-i))oo{a;2(l+'-)z; a;*, x2('-i))o 
(aJ5; 5 X )oo ' 

(a^2r-l^ 3;2r)^ 

(a;2; x'^'^)oo ' 
(a;-i2;a;2{'-i))oo ' 



(A.2) 



B Proof of (4.6) and (4.7) for e = — and general k 



In this appendix we shall give an inductive proof of ( |4.6| ) and (4.7) for e = — and general k. The 
first step is to show that if we have 

k 



u]^±{0\k,k)% = ^±{C-^)\k,k)% 



for a given value of /c, then it follows that 

^-{0% k + l)% = ^-{C')\k, k + 1)% 



< V k + l 



(B.l) 



(B.2) 



is also true for the same value of k. Using ( 4.29| ) and the vertex operator commutation relations 
( |3.49D , we can rewrite ( |B.2| ) as 



u] < W 



k 



k + l 



k + l k 



C+l-u] $_(x2^+2)$+(C) 



+w 

{"'I 

+W 



k k-1 

k + l k 

k k + l 

k + l k 

k k-1 

k + l k 



c+l 
B 



c+l-uj «>+(2;2^+2)$_(C)| \k,k)% 

C+l+u] «>_(x2^+2)^.+ (C-l) 



C+l + uU+(x2^+2)$_(C-l) \k,k)''^\ 



(B.3) 



We can now use the identities 



vk+1 




V k + l 





k k + l 

k + l k 

k k-1 

k + l k 



c + 1 — u 



c + 1 — u 



( A: + 1 



K 



(c+l) 



k-1 



u] W 



k 



k + l 



u] W 



k+l k 

k k-1 
k+l k 



c+l + u 



c + l + u 



and the given relations ( |B.1| ) to re-write the LHS of ( |B.3| ) such that it is identically equal to the 
RHS. 
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The next step in the proof is show that if 



(^A: - 1^1 cI,_(C)|fc _ 1, k)% = $_(r')|A: - 1, k) 



B 



(B.4) 



is true for a given k, then 



K 



k-1 



u]<i>^{0\k,k)% = ^^{C-')\k,k)% 



(B.5) 



follows. To show this we rewrite (B.5), using ( 4.30| ) and the commutation relations ( 3.54| ), as 



Res^^_,^+2cx{aC)K^> {k -ll\u) ^;(e)1>-(C)|A: - l,k) 

= Res^=_,i+2cx(cc)n(0^-(r')i^ - 1, 



B 



(B.6) 



The equality of the two sides is then established using (B.5) and the identity 



(B.7) 



We have now completed the inductive steps. In Section 4 we proved (^|^) for e = — and k = \ 
(( [4.6[ ) with e = + and general k is true by construction). The proof of ( [4.6[ ) and (4.7) is complete. 
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